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ABSTRACT 

The abundances of metals in the intergalactic medium (IGM) can be used to constrain the amount of star 
formation at high redshift and the spectral shape of the ionizing background radiation. For both purposes it is 
essential to measure the abundances in regions of low density, away from local sources of metals and ionizing 
photons. Here we report the first detection of O VI in the low-density IGM at high redshift. We perform a pixel- 
by-pixel search for O VI absorption in eight high quality quasar spectra spanning the redshift range z = 2.0-4.5. 
At 2 <J z 3, we clearly detect O VI in the form of a positive correlation between the H I Lya optical depth and 
the optical depth in the corresponding O VI pixel, down to thi ^ 10"'. This is an order of magnitude lower in thi 
than the best C IV measurements can probe and constitutes the first clear detection of metals in underdense gas. 
The non-detection of O VI at z > 3 is consistent with the enhanced photoionization from a hardening of the UV 
background below z ^ 3 but could also be caused by the high level of contamination from Ly series lines. 

Subject headings: cosmology: observations — galaxies: formation — intergalactic medium — quasars: 
absorption lines 



1. INTRODUCTION 

The abundance of metals in the low-density intergalactic 
medium (IGM) at high redshift is of interest for at least two rea- 
sons. Firstly, the metallicity can be used to distinguish between 
different enrichment mechanisms and to constrain the amount 
of star formation at high redshift. Secondly, the relative abun- 
dances of ions with different ionization potentials can be used 
to derive the spectrum of the integrated UV background from 
stars and quasars. For both purposes it is essential to measure 
the metal abundances in gas of low density, away from the influ- 
ence of local sources of ionizing radiation and locally produced 
metals. 

Resonant Lya absorption by neutral hydrogen along the line 
of sight to distant quasars results in a forest of absorption lines 
bluewards of the quasar's Lya emission line. With the advent 
of the High-Resolution Echelle Spectrograph (HIRES) on the 
Keck telescope, it became clear that many of the high column 
density Lya lines (A^hi ^ 10'^'^ cm"^) show associated C IV ab- 
sorption (e.g., Cowie et al. 1995). Although increasing the 
signal-to-noise ratio (S/N) reveals C IV systems of progres- 
sively lower column density, C IV lines corresponding to low 
column density Lya lines are generally far too weak to detect 
(e.g., EUison et al. 2000). 

The success of cosmological simulations in reproducing the 
observations of the Lya forest has convincingly shown that the 
low column density forest (A^hi < 10'"* ''cm"^) at high redshift 
(z ^ 2) arises in a smoothly fluctuating IGM, with individual 
lines corresponding to local density maxima of moderate over- 
density ip/p <, 10). The interplay between photoionization 
heating and adiabatic cooling that is due to the universal ex- 
pansion results in a tight temperature-density relation, which 
is well approximated by a power-law, T = To(p/p)''~^ (Hui & 
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Gnedin 1997). The Lya optical depth is proportional to the 
neutral hydrogen density, which in photoionization equilibrium 
is proportional to p^T""^^. Hence the optical depth depends on 
the underlying baryon overdensity, thi oc (p/ pf-^^-''^'' . The 
constant of proportionality is an increasing function of redshift, 
at z ~ 3, an optical depth of one corresponds to slightly over- 
dense gas. 

Cowie & Songaila (1998) realized that the tight correlation 
between optical depth and gas overdensity makes a pixel-by- 
pixel analysis of metal abundances useful. By measuring the 
median C IV/H I optical depth ratio as a function of thi, they 
were able to show that the IGM is enriched with a roughly con- 
stant C IV/H I ratio down to thi ~ 1 at z 3. Ellison et al. (2000) 
showed explicitly that directly detected C IV lines account for 
only a small fraction of the total amount of metals in the forest, 
as inferred from a pixel analysis. 

Although other metal lines have been detected, C IV 
(AA1548, 1551) has so far proven the most sensitive since it 
is the strongest line redwards of Lya (A1216), where there is 
little contamination from other absorption lines. However, sim- 
ulations show that at high redshift, photoionized O VI (AA1032, 
1038) should be a much better probe of the metallicity of the 
IGM in regions close to the mean density (Chaffee et al. 1986; 
Rauch, Haehnelt, & Steinmetz 1997; Hellsten et al. 1998). In 
practice it has proven very difficult to detect O VI at high red- 
shift because it lies deep in the Lya and Ly/3 (A1026) forest. 
Since Lu & Savage (1993) established the presence of O VI 
in intervening absorbers by stacking Lyman limit systems that 
show C IV absorption, O VI has been detected at high redshift 
in several Lyman limit systems (Vogel & Reimers 1995, Kirk- 
man & Tytler 1997; 1999). Dave et al. (1998) (see also Cowie 
& Songaila 1998) carried out a thorough search for O VI but 
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found only evidence for O VI associated with C IV absorption 
in high density gas. 

At low redshift (z ^ 1), O VI absorbers have been shown to 
be common (e.g., Bergeron et al. 1994), so common in fact 
that they may harbor a large fraction of the baryons (Buries &, 
Tytler 1996; Tripp, Savage & Jenkins 2000). These low red- 
shift systems probably correspond to hot, shock heated gas in 
the potential wells of (groups of) galaxies, which may well be 
coUisionally ionized. 

We have carried out a pixel-by-pixel search for O VI in the 
spectra of eight quasars, spanning the redshift range z ~ 2.0- 
4.5. We clearly detect O VI at 2 ^ z ^ 3, down to very low 
H I optical depths, thi ^ 10"', which is an order of magnitude 
lower than the best C IV measurements can probe. 



TABLE 1 

List OF samples 



Sample 


Zmin 


Zmax 


(S/N)ovi 


(S/N)hi 


QSO 


Zem 


1101 


1.97 


2.10 


21 


78 


Ql 101-264 


2.14 


2233 


2.07 


2.18 


12 


30 


J2233-606 


2.24 


1122 


2.03 


2.34 


27 


66 


Q1122-165 


2.40 


1442 


2.51 


2.63 


8 


61 


Q1442+293 


2.67 


1107 


2.71 


2.95 


10 


49 


Q 1107+485 


3.00 


1425 


2.92 


3.14 


42 


113 


Q 1425+604 


3.20 


1422 


3.22 


3.53 


46 


105 


Q1422+231 


3.62 


2237 


4.15 


4.43 


31 


32 


Q2237-061 


4.55 



2. OBSERVATIONS AND SAMPLE DEFINITION 

We analyzed spectra of the eight quasars listed in column 6 
of Table 1 . The spectra of Q 1 1 1 , J223 3 and Q 1 1 22 were taken 
during the Commissioning I and Science Verification observa- 
tions of the UV- Visual Echelle Spectrograph (UVES) and have 
been released by ESO for public use. They were reduced with 
the ESO-maintained MIDAS ECHELLE package (see Kim, 
D'Odorico, & Cristiani 2000, in preparation) for details on 
the data reduction). The other spectra were obtained with the 
HIRES spectrograph (Vogt et al. 1994) on the Keck telescope. 
The reduction procedures for these quasars can be found in Bar- 
low & Sargent (1997) and Sargent, Barlow & Rauch (2000, in 
preparation). The HIRES and UVES spectra have a nominal ve- 
locity resolution of 6.6 and 6.7 kms"^ (FWHM) and a pixel size 
of 0.04 and 0.05 A respectively. To avoid confusion with the 
Ly/3 forest, we only used the Lya region between the quasar's 
Lya and Ly/3 emission lines. In addition, spectral regions close 
to the quasar (typically 5000 kms~') were excluded to avoid 
proximity effects. Regions in the Lya forest thought to be con- 
taminated by metal lines were excluded. 

The mean absorption increases rapidly bluewards of the 
quasar's Ly/3 emission Une becaue of the growing number of 
absorption lines other than Lya. For quasars at z 3 the in- 
crease in the absorption is dominated by higher order Ly lines, 
corresponding to Lya lines which fall in between the quasar's 
Ly/3 and Lya emission lines (i.e. in the Lya forest). At lower 
redshift the hydrogen lines have smaller column densities and 
metal lines contribute significantly to the mean absorption. The 
increase in the mean absorption towards the blue, raises the 
threshold for detecting O VI absorption in the Lya forest con- 
siderably. Moreover, the S/N in the O VI region generally de- 
creases towards the blue, raising the detection threshold further. 
We therefore analyzed the red and blue halves, i.e. the high 
and low redshift halves, of each Lya forest spectrum separately 
(for those quasars for which our spectral coverage extends far 
enough into the blue). This turned out to be crucial, as we were 
unable to detect O VI in any of the low redshift half samples. 
All the high redshift half samples are listed in Table 1 . 

*The spectrum is divided into chunks of 2 A and 250 spectra of length equal to the original spectrum are generated by picking chunks at random (with replacement, 
each chunk can be picked more than once). The la error in the median is then the square root of the variance in the median over the bootstrap realizations. 



3. METHOD 

We search for O VI using a variant of the pixel technique in- 
troduced by Cowie & Songaila (1998). We measure the optical 
depth in each Lya forest pixel and in the corresponding O VI 
pixel. The H I-O VI pixel pairs are binned according to their 
H I optical depth and the median H I and O VI optical depths 
are computed for each bin. This technique enables us to probe 
the O VI abundance down to much lower densities than is pos- 
sible by fitting lines. Its main advantage over, for example, the 
stacking method is its robustness. Being a median, it is rela- 
tively insensitive to (non-Gaussian) noise, contamination from 
non-0 VI lines, problems that severely compromise the stacking 
method, which is a mean (see Elhson et al. 2000 for a critical as- 
sessment of both the stacking technique and the pixel technique 
for the case of C IV). Its main limitation is that it measures an 
apparent O VI optical depth, Tovi.app, which is in general not the 
same as the true O VI optical depth because of contamination. 

The effect of contaminating fines is reduced by defining 
TovLapp = min(Tovi, 1032, 2tovi, 1038), where we used the fact that 
(/-^)i032 = 2(/A)io38. Taking the minimum doublet component 
does of course not guarantee that Tovi.app = tovi. For some pix- 
els both components may be contaminated, which would re- 
sult in an overestimate of the O VI optical depth. However, if 
the contamination is negligible, then noise will cause us to un- 
derestimate Tovi- For C IV, which falls redwards of the Lya 
forest, noise is the limiting factor and Ellison et al. (2000) 
therefore only used the weaker component if a 3<j detection 
was predicted, based on the signal in the stronger compo- 
nent. For O VI however, contamination is a much greater prob- 
lem and we therefore take the the minimum if the expected 
signal in the weaker component is greater than 0.5cr, i.e. if 
exp(-0.5rovi, 1032) < 1 -0.5crovi, 1038, where ctovi, io38 is the rms 
amplitude of the noise at the position of O VI, 1038 (as given by 
the normalized noise array). 

When a pixel is close to saturation [exp(-rHi) < 0.5ct], we 
use up to 10 higher order Ly lines to determine the optical 
depth: TLya = min(rLy„/LyQALya//Ly«ALy„). We limit the effect 
of noise features by requiring 0.5(t„ < exp(-TLy„) < 1 -0.5(t„ 
(i7„ is the rms noise amplitude at the position of Lyn). The num- 
ber of higher order lines available is different for each quasar 
spectrum and varies also within a single spectrum. However, 
since the number of pixels with a given rm decreases rapidly 
for Thi > 1, the misbinning of H I due to missing higher orders 
affects only the highest thi bins. 

We use logarithmic thi bins of size 0.35 dex (0.5 dex for the 
two lowest redshift quasars) and estimate the error in the medi- 
ans by bootstrap resampUng^. 
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Fig. 1 . — Results from the pixel analysis. For all HI hya pixels within the redshift range indicated in each panel, the corresponding OVI pixel was determined 
and the optical depths were measured. Pixels corresponding to higher order Ly lines and to the weaker line of the OVI doublet were used as described in section W. 
HI-OVI pixel pairs were binned according to their HI optical depth and the data points indicate the median HI and OVI optical depth for each bin. The measurea, 
apparent OVI optical depths include contributions from noise and contaminating non-OVI lines. The dashed lines show the median, apparent OVI optical depth for 
the full set of HI pixels. If the contribution of OVI to the median absorption is negligible, then the dashed lines are effectively the OVI detection limits. Vertical error 
bars are Icr en'ors, horizontal error bars are smaller than the symbols and are not shown. For z < 3 rovi.app and rm are clearly correlated, down to optical depths as 

low as THi ~ 10"' . A correlation between Tovi.app ™d 'm implies that OVI absorption has been detected in the Lya forest. 



4. RESULTS 

Figure 1 shows the results from a pixel analysis of the sam- 
ples listed in Table 1 . The apparent O VI optical depth is clearly 
correlated with thi at z < 3 down to thi ^ 10~\ indicating that 
O VI is detected in the Lya forest. At z ~ 3.2-3.5 we detect O VI 
only at very high H I optical depths, thi > 10, while at z > 4 rovi 
and Thi appear uncorrelated. The transition seems to occur at 
z ^ 3.0-3.5. If contamination and noise, which are independent 
of Thi, dominate over absorption from O VI, then Tovi.app will 
flatten off at the detection limit. Ellison et al. (2000) analyzed 
a very high quality spectrum of Q1422+231 ((S/N)civ 200) 
and found that rciv flattens off below thi ~ 2. We find that the 
correlation between Tovi.app and thi continues down to much 
lower Thi, indicating that oxygen has been detected in gas of 
very low density contrast (p/p < 1). 

The horizontal, dashed lines indicate the median Tovi.app cor- 
responding to a random H I pixel, i.e. the median Tovi.app for the 
set of all H I pixels. This reference level includes contributions 
from O VI, from contaminating non-O VI lines and from noise. 
For z ^ 3, the contribution of O VI to this level is negligible 



and the level is effectively the detection limit. For z ^ 3 on the 
other hand, O VI accounts for a significant part of the total ab- 
sorption. Hence the detection limit is lower than the reference 
level and Tovi.app drops below the dashed line at low thi. 

We stress that because of noise and especially contamination, 
the measured, apparent O VI optical depths could differ consid- 
erably from the true values. We tried to estimate the effects of 
noise and contamination, which vary systematically across the 
spectrum, by analyzing synthetic spectra generated by drawing 
H I Voigt profiles at random from observed line lists until the 
mean absorption in the Lya forest equals the observed value. 
We then added the corresponding O VI absorption lines using 
a fixed O VI/H I column density ratio and a fixed /j-parameter 
ratio. Finally, higher order Ly lines and noise were added using 
the H I line list and the noise array respectively. We found that 
the simulations agree reasonably well with the observations for 
A^(0 VI)/A^(H I) 10"^-10"', although the results ai-e sensitive 
to the assumed ratio of /^-parameters. Although encouraging, 
these Monte Carlo simulations are clearly too simplistic to de- 
rive the intrinsic O VI/H I ratios. The effects of noise and con- 
tamination can only be modeled accurately by calibrating the 
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method against synthetic spectra extracted from hydrodynamic 
simulations. This should be done separately for each observed 
spectrum, using synthetic spectra that resemble the observed 
spectrum in detail. We leave such a quantitative analysis for a 
future paper but note that the effects that prevent us from de- 
riving the actual O VI abundance, will reduce any intrinsic cor- 
relation between thi and rovi,app. making a detection of such a 
correlation a robust result. 

For a fixed O VI abundance, the value of the reference level 
(Fig. 1, dashed lines) increases with decreasing S/N and with 
increasing redshift. The latter effect is more important and is 
due to the increase in the number and strength of contaminat- 
ing Ly series Unes with redshift, which is a direct consequence 
of the expansion of the universe. As the universe expands, the 
column densities of lines corresponding to absorbers of a fixed 
overdensity decrease (the evolution of the ionizing background 
and the growth of structure are less important). This results in 
an increase of the amount of contamination, and thus the detec- 
tion threshold, with redshift. Note that this effect is opposite 
to the increase in contamination towards the blue (i.e. low red- 
shift) end of a single quasar spectrum, which is a consequence 
of the growing number of interlopers (see section 2). Like the 
column density, the H I optical depth corresponding to a fixed 
overdensity also increases with redshift. This complicates the 
comparison of the thi level below which the correlation with 
O VI flattens off in samples of different redshifts. While an H I 
optical depth of one corresponds to about the mean density at 
z ^ 3, it corresponds to an overdensity of a few at z ^ 2 and to 
slightly underdense gas at z ~ 4. 

For the quasars Ql 107, Q1422 and Q2233 our spectral cov- 
erage is sufficient to analyze the blue half of the forest. We find 
no correlation between tovi and thi in any of these samples. 
As discussed above, this is a consequence of the large number 
of contaminating absorption lines from intervening absorbers 
and/or the low S/N of the data and does not imply that there is 
no O VI. In principle, even more O VI could be uncovered by 
looking at a smaller part of the spectrum. The contamination 
in the upper redshift quarter for example, will be lower than in 
the upper half. However, decreasing the number of pixels fur- 
ther would bring the redshift path into the regime of large-scale 
structures and the samples would therefore no longer be rep- 
resentative. This may already be the case for the two lowest 
redshift samples, for which the number of pixels showing H I 
absorption is much smaller than for the higher redshift samples. 
Nevertheless, we have tried smaller samples and found that the 
correlation strengthens somewhat for 1422 but not for 2237. 

Figure 2 demonstrates that the O VI that we detect at low 
Thi is not associated with detectable C IV absorption. The solid 
line indicates the apparent O VI absorption in 1442. The dashed 
line shows the result obtained after all pixels that have nonzero 
C IV absorption at the 1 g level [exp(-Tciv) < 1 - Icr] have been 
excluded. The dashed curve only falls below the solid curve for 
Thi ^ 2, indicating that some of the stronger O VI systems have 
been excluded. The fact that the two curves are almost indistin- 
guishable for Thi ^ 2 shows that O VI is detectable in absorbers 
for which C IV absorption is too weak to detect and is thus 
direct observational proof of the prediction that O VI is a more 
sensitive probe of the metaUicity in low-density gas than C IV. 



1442 




median(XLya) 

Fig. 2. — The effect of excluding CIV absorbers. The solid curve is for sam- 
ple 1442. The dashed curve shows the result of the pixel analysis after having 
excluded all pixels for which the corresponding CIV absorption is detectable. 
The curves only start to differ at thi > 2, indicating that only the strong OVI 
absorbers show significant CIV absorption. 

5. DISCUSSION 

We have searched for O VI absorption in the spectra of eight 
quasars spanning the redshift range z ^ 2.0-4.5 and reported 
the first detection of O VI in the low-density IGM at high red- 
shift. We analyzed the spectra pixel by pixel and detected a 
positive correlation between the optical depths in H I and in the 
corresponding O VI pixels down to thi ~ 10"' at 2 < z < 3. 
This is an order of magnitude lower in thi than the best C IV 
measurements can probe and constitutes the first firm detection 
of metals in underdense gas. We showed that the O VI signal 
does not come from systems for which C IV is detectable, and 
this confirms the prediction that O VI is a much better tracer of 
metals in low-density gas than C IV. 

Although the very strong absorbers may well be coUisionally 
ionized (e.g., Kirkman & Tytler 1999), the volume filling fac- 
tor of this hot gas is much smaller than that of photoionized gas 
(e.g., Cen & Ostriker 1999). The ubiquitous O VI detected here, 
arises in gas that is too dilute for collisional ionization to be ef- 
fective. Photons of at least 8.4 Ryd are needed to create O VI 
and hence we would not expect to detect photoionized O VI if 
the reionization of He II, for which photons of at least 4 Ryd 
are required, is not yet complete. Measurements of the He II 
Lya opacity (Heap et al. 2000 and references therein) and the 
temperature of the IGM (Schaye et al. 2000, Ricotti, Gnedin & 
ShuU 2000) suggest that the reionization of helium may have 
been completed at z ^ 3. Although our non-detection of O VI at 
z > 3 is consistent with enhanced photoionization from a hard- 
ening of the ionizing background below z 3, it could also be 
caused by the increased level of contamination from Ly series 
lines. 

Our detection of O VI down to thi ~ 10 ' shows that the IGM 
is enriched down to much lower overdensities than have been 
probed up till now. A comparison with simulations could clar- 
ify how noise and contamination from Ly series lines affect the 
apparent O VI optical depth, what fraction of the low-density 
gas is enriched and what the scatter in the O VI/H I ratio is, 
thereby providing strong constraints on theoretical models for 
the enrichment of the IGM. 
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